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Abstract. We derive asymptotic expansions for the displacement at the bound- 
ary of a smooth, elastic body in the presence of small inhomogeneities. Both 
the body and the inclusions are allowed to be anisotropic. This work ex- 
tends prior work of CapdeBoscq and Vogelius (Math. Modelling Num. Anal. 
37, 2003) for the conductivity case. In particular, we obtain an asymptotic 
expansion of the difference between the displacements at the boundary with 
and without inclusions, under Neumann boundary conditions, to first order 
in the measure of the inclusions. We impose no geometric conditions on the 
inclusions, which need only be measurable sets. The first-order correction con- 
tains an elastic moment tensor M that encodes the effect of the inclusions. In 
the case of thin, strip-like, planar inhomogeneities we obtain a formula for 
M only in terms of the elasticity tensors, which we assume strongly convex, 
their inverses, and a frame on the curve that supports the inclusion. We prove 
uniqueness of M in this setting and recover the formula previously obtained 
by Beretta and Prancini (SIAM J. Math. Anal, 38, 2006). 



1. Introduction 

Let f2 C M. d , d > 2 be a smooth bounded domain representing the region occupied 
by an elastic body. Let Co = Co (a;) be a smooth background elasticity tensor in O. 
Let w e C O be a set of measurable small inhomogeneities and let Ci = Ci(x) be the 
smooth elasticity tensor inside the inhomogeneities. Let ip g H^ 1 ^ 2 {dVl) represents 
a traction on dVt and U the corresponding background displacement field which 
satisfies the system of linearized elasticity: 

div(CoW) = in fi 

(c v[/> = v on an, 

Let 

C £ = CoXn\ w< +Cix We , 
and consider the perturbed displacement field solution to 

div(C £ Vu £ ) = in 
(C e Vw e )^ = if on dn. 

One goal of this paper is to obtain an asymptotic formula for u e — U on the boundary 
of Q as the measure of lo c approaches zero. The formula we derive generalizes 
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those already available in case of homogeneous isotropic bodies with diametrically 
small (see [3], and, also p], [2]) or for thin (see [4]) inhomogeneities. To derive 
the asymptotic expansion we follow the approach introduced by Capdeboscq and 
Vogelius in j6] (see also [7]) for the conductivity equation and we establish a formula 
in the case of arbitrary elastic tensors Co and Ci. More precisely we show that for 
yedn 

(1.1) (u en - U)(y) = K | f M(x)VU(x) : V7V(x, y)d l i x + o(|cj en |). 

along a sequence of configurations {w £ „} whose measure tends to 0. Here, N is the 
Neumann function corresponding to the operator dw(CoV-), fj, is a Radon mea- 
sure, the elastic moment tensor M 6 L 2 (il, dfi) and VU represents the symmetric 
deformation tensor. 

For particular geometries like diametrically small or thin inhomogeneities the as- 
ymptotic expansion holds for (u e — U)(y) as e — > and one can characterize the 
measure /i and the tensor M. In particular, if w e = z + eB, where the center z G f2 
and B is a bounded domain, then fi is a Dirac function concentrated at z. If fur- 
ther both Co and Ci are homogeneous and isotropic, the tensor M can be explicitly 
computed and carries information about the geometry of B and about the elastic 
parameters of Co and Ci ([3]). If uj £ — {x G il, dist(x, cr ) < e}, where a is a 
simple smooth open curve in the plane, [i reduces to a Dirac measure supported 
on <7o- If again the phases are isotropic, M can be explicitly determined by the 
transmission conditions for u e [3]. 

In the second part of the paper we analyze the case of thin inhomogeneities in a 
planar domain in the case of arbitrary elasticity tensors. 

In the case of isotropic homogeneous tensors, the idea, used in [4] to derive the 
asymptotic expansion for u e — U , is to apply fine regularity results for solutions 
of elliptic systems with discontinuous coefficients by Y.Y. Li and L. Nirenberg [TU] 
and to use the transmission conditions to derive the tensor M which satisfies 

(1.2) (Ci - Co)Vt**(a;) = M(x)Wu e t (x), 

whereas Vu* and Vwf denote the values of the deformation tensor inside and outside 
the inclusion at a point x on its boundary. Note that the deformation tensors are 
related by transmission conditions across duj e . 

One of the difficulties we encountered in deriving the expansion in the anisotropic 
case is the the direct derivation of fj 1 . 2|) from the transmission conditions. To 
construct M, we follow the work of Francfort and Murat [8] on the calculation of 
the effective properties of laminated 2-phase elastic composites. Indeed, one can 
view moment tensors as limits of effective tensors as the volume fraction of one of 
the phases tends to 0. 

The paper is organized as follows. In section 2, we state a general representation 
formula of the form (jl.ip for anisotropic elastic inhomogeneities embedded in an 
anisotropic background medium. The asymptotic expansion is proved in section 3. 
Properties of the elastic moment tensor M are established in section 4. In section 5, 
L0 e is assumed to be a thin strip-like planar inclusion. Firstly, relying on the uniform 
Holder regularity of u £ and on Meyer's theorem, we give a direct derivation of the 
asymptotic expansion similar to that in [4] , under the assumption that there exists 
a tensor M, independent of e, that satisfies (|1.2|) . Secondly, we prove existence of 
such M, invoking the result of Francfort and Murat mentioned above [8! . Thirdly, 
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we show that the asymptotic expansion of theorem 2.1 coincides with that obtained 
in theorem l5.ll Finally, in the appendix, we recall classical regularity results for the 
system of elasticity, and we prove how Caccioppoli inequality and Meyer's theorem 
also hold for the system of elasticity. 
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2. Notations, assumptions and main result 

Let f2 C M. d , d > 2 be a bounded, smooth domain. For x £ <90, let us denote by 
v(x) the normal direction to dQ at point x. We use the following notation: 

Notation. Let C be a A-th order tensor, let A and B be d x d matrices, and let u, 
v denote vectors in M. d . We set: 

u-v = Yf j =iUjV j Ac V'' , A,,r, €.A = Yf k ,i=i <C m A H 

(CA) v = E^,;=i CijkiAkiVj CA:B = J2ij.k.i=i ^ijkiAkiBij, 

l^l = (E«4-) 1/2 - 

Moreover, we denote by A = (A + A T )/2 the symmetrization of the matrix A. In 
particular, given a vector valued function u defined in fl, we denote by Vu the strain 
Vu = | (Vu + (Vu) r j . 

Let Co £ C 1:a (ft), for some a £ (0,1), be a fourth order elasticity tensor that 
satisfies the full symmetry properties: 

(2.1) (Ca(x))ijki = (Ca(x))kUj = (Co(^))j»fci VI < k,l < d and ie!l, 
and the strong convexity condition, i.e., there exists a constant Ao > such that 

(2.2) Cq(x)A : A > Ao|A| 2 for every d x d symmetric matrix A and x £ Q. 
Let i}j £ H^ 1 / 2 (dfl) satisfying the compatibility condition 



(2.3) / ip-R = 0, 

Jan 

for every infinitesimal rigid motion R, that is R{x) = Wx + c for some skew- 
symmetric matrix W and c £ M. d . 

The background displacement field U £ H(fl) is defined as the solution to 

f div(CoW) = Oinr! 
{ ' ' 1 (C Q VU)iy = V on dtt, 

where H(Cl) is the space of vector valued functions given by 

H (fi) = \u £ H 1 (n;M d ) such that / u da = 0, / (Vu - (Vu) T ) dx = 1 
I Jan Jn J 
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Let uj e denote the a subset of f2, that contains one or several inhomogeneities. We 
assume that w e is measurable and separated from the boundary, that is d(u> e , dtt) > 
do > 0. We also assume that the measure \u e \ > tends to as e 0. 
Let Ci denote the elasticity tensor inside to e . We assume that Ci € C a (Vt) is fully 
symmetric and strongly convex, i.e. 

(2.5) C\{x)A : A > Ao|A| 2 , for every d x d symmetric matrix A and x G ft. 
Let C e be the elasticity tensor in the presence of the inhomogeneity 

(2.6) C e = CoXn\ We +Cix We , 

and consider the corresponding displacement field u t € H(Q) solution to 



(2.7) 



div(C e Vw e ) = in 
(C e Vu £ )^ = ip on cTJ . 



For existence and uniqueness of solutions to (|2.4p and (|2.7[) in H(£l) we refer to 
pj], for example. 

Since d(w £ ,<9f2) > do > 0, there exists a compact set Kq, independent of e, such 
that 

(2.8) w £ C K Q C n and dist(w £ , Q \ K ) > d /2 > 0. 

We also introduce the Neumann matrix for the operator div(CoV-), i.e. the weak 
solution to 



f div(C V7V(-,y)) = -6 y I d mn 
\ (C VN(-,y)u = -^^ondn, 

that satisfies the normalization conditions 



(2.9) 



(2.10) / N{x,y)da x = 0, / {V x N(x, y) - V x N{x, y) 1 ) dx = 0, 
Jan Jn 

where 1^ is the d-dimensional identity matrix. 

For the existence of such Neumann matrix and its behavior for x close to y we 
refer to [9] where existence and regularity of the Green's matrix for weakly elliptic 
systems is considered. 

The following result generalizes the compactness result of [5] to the case of elastic 
inclusions: 

Theorem 2.1. Let uj en be a sequence of measurable subsets satisfying \2. 8\) such 
that, as n — > oo ; |w £ J — > and 

(2.11) \oj €n \~ 1 Xu>i n dx — > d/i in the weak* topology o/(C(fi))', 

for some regular positive Borel measure fi, such that J^dfi = 1. 
Given i)j € H~ 1 / 2 (dQ) satisfying \2. 3\) . let U and u tn denote the solutions to \2.J$ 
and {2. 7| ) respectively. There exists a subsequence, not relabeled, and a fourth order 
tensor M £ L 2 (tt, d/i) such that, for y e dfl, 



(2.12) (u eB -Z7)(i/) = |w eB | / M(x)VU(x):VN(x,y)dfi x + o(\oj 6n \). 

We will prove this result in the next section. 
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3. Proof of theorem 12.11 

3.1. Preliminary estimates. Let F G and / G H^ 1 ^ 2 ^) satisfying the 

compatibility conditions 

F dx = f f da x and / F ■ Rdx = / f ■ Rda x , 
Jan Jo. Jan 

for every infinitesimal rigid motion R. 

Let V and v e in H(Q) solve 



(3.1) 

and 

(3.2) 

respectively. 



div(C W) = F in fi 

(<c vv> = f on an, 



div(C e V'u £ ) = F in Q 
(CqVuc)^ = /on dQ, 



Lemma 3.1. Let F e C a (0), im'tt < a < 1 and let < n < 1/d. There exists a 
constant C > 0, smc/i i/iai 

(3.3) |K - V\\ m{n) < C\u e \ 1/2 (||F|| c -(n) + \\F\\ 
and 

(3.4) ||t, e - V\\ L * m < Cl^+i-v (||F|| ca(n) + ||F|| ff -i (n) + H/Hff-x/^ao)) ■ 

Proof. We adapt the arguments of ([6]) to the system of elasticity. Since V and v e 
solve (|3.1|) and (|3.2I) respectively, for every iu g if 1 (fi), we have 

/ C e V(w £ - y) : Vffl dx = / (Co - C X )W : Vw dx. 
By choosing w = v e — V and applying Korn's inequality, we show that 

|v(« e - v)\ 2 dx < c [ |v(v e - v)\ 2 dx < cm 1 ' 2 !! w|U- K) || v(« £ - v)|U a(n) . 



It follows from interior regularity results for the elasticity system with regular co- 
efficients (see, for example, Theorem 6. Ill, chapter 2 in [5]) that 

l|W|Uoo (we) < C (\\V\\ HHQ) + ||F|| C a (Jfo) ) 

(3-5) < C i\\F\\ H - 1(fl) + H/Hff-v^Q) + ||F||c«(K )) , 

and, hence, 

We also have, for any w E i? 1 (fi), 

(3.6) { C V(w e - V) : \7wdx = f (C - Ci)Vw e : Vw dx. 

Let us select w £ H(fl) as the solution to 

j div(C Vw) = V - v e in Q, 

1 (CqVw)^ = J n (V — v e )dx on <9fi, 
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By the smoothness assumption on Co and by interior regularity estimates (see, for 
example, Theorem 2.1, chapter 2 in [5]) we have that 

IM|iP(if D ) < C{\\v t - V\\ L 2 {n) + \\w\\ H i {n) ). 

By Korn and Poincare inequalities 

IM|jri(n) < 0\\Vw\\ L3(n) <C||VHU»(ii) < \\V - v e \\ L 2 {n) 

and, hence, 

(3-7) \\w\\ HHKo) <C\\v € -V\\ L 2 m 

By Sobolev Embedding Theorem, we have that Vic € L p (Kq) for every 1 < p < d* 
where d* = -j^ f° r d > 2 and d* = +oo for d — 2, and 

(3.8) Qf_ \Vw\*dx^ P < Cp|Mlfl*(jr ) < C\\v e - V\\ L n n y 

Let us choose q e (^pj , 2) and p such that | + | = 1. Notice that, p G (1, gP). By 
inserting iu into (|3.6[) we obtain 

(v e -V) 2 dx = / C V(w e - V) : Vwdx = I (C - C^Vv, : Vw dx 



(3.9) < C \Vv e \i 



(3.10) < C / |Vu e |« ||w e - ^IU» ( n) 



Now, by Holder inequality and (|3.5[) we get 

||V» e |U« (ll) .) < ||V(U £ - V)|U« (u;e) + ||W||i, (u .) 

(3-11) < |w e |i-*||V(w e - V)]!^^.) + |w e |«||VV]Uoo (Wjl) 

< C^l 1 /" (||F|U-i (n) + ||/|| ff -i/ a(8 n) + II^IU-pi)) 

A combination of (f3~9|) . ff378]l and ([3~TT|) yields 

K - V|| £ »(n) < C q \u) e \ 1 l q (||F|| 

Note that since for q \ we have - /• ^ + |. It follows that given any 

< 77 < 1/d there exists a constant C such that 

Ik - ^IU>(n) < C|« e |i + 4-" (I^IU-Mn) + ll/llH-v 2(9 n) + ll^llc=(o)) 

□ 

3.2. Definition of the elastic moment tensor. Let 

(3.12) v %3 = i (eiXj + QjXi) - Cij 

where e^ is the i-th coordinate direction and 

ClJ = 2\dn\ j an ( - GlXj + GjXi ^ d<J 
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and consider v l J £ H(£l) solution to the problem 
(3.13) 



div(C e Vwf ) 
(C e Vt)* J > 



div(C Vz/ J ') in Q 
(C W J > on 90, 



Observe now that 



i—rX-.tCi- Co) V«« 

We 



< 



L!(n) 



1 



C )V(«« - « <3 ') 



da; 



- Co)V«*- 



dx. 



By Lemma 13.11 and recalling that Co and Ci are bounded and that 



(3.14) 
we have 



Vw lJ = Vv ij = - fa <E) e 3 + e 3 ® e*), 



T-rWCi-QOVi;? 



< C. 



L x (n) 



Hence, possibly extracting a subsequence, we may assume that 
(3-15) k„|x Wen (Ci - C )Vw« -> <M ijlm , 

in the weak* topology of C°(f2), where dMijj m is a regular Borel Measures with 
support in Kq. Let $ G C°(f2). By dehnition of dMy; TO , we see that 



$dM 



< lim„ 



- lim - — 



ijlm 



1 



lim 



\uj r _ 



(Ci - Co)V(«* j - f ij ') |*| dx 



(Ci - C Q )Vv y 



< linu 



Kl 1/a 

1/2 



|$| dx 

1/2 



|V(v£ -v ij )\ 2 dx 

1/2 



O Pen I 



rX^J$| 2 dx 



1/2 



< c 



\*\ J dn 



Hence 



$dM 



ijlm 



is a bounded functional on L 2 (f2, d/x), and 



QdMijim = / $Mi,- im d/i 



for some function Mjj/ m € L 2 (f2,d/x). The tensor M actually relates to the weak 
limit of u e , as the next lemma expresses: 

Lemma 3.2. Let U andu e denote the solutions to \2.1$ and \2. 7\ ) for ip g i7 _1 / 2 (<9f2) 
satisfying the compatibility conditions \2. 3\) . Let uj erl such that |w Cn | — > be a se- 
quence for which S2.8\) . \2.11}) and S3.15\) hold. 
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Then, ^ ; Xu> Cn (Ci — <Co)Vu Cn dx is convergent in the weak* topology of (C°(f2))' 
with 

lim r^-rXw, (Ci - Co)Vw e „cfe = MWd/i 

Proof. It suffices to prove that we may extract a subsequence of {w e „} such that 
r-p r Xu„ (Ci — Co)Vu e _.dx converges to MVUdu. The fact that the limit is 

independent of the particular subsequence guarantees that the entire sequence is 
convergent. 

Proceeding as for v l J , we see that 



1 



-X^ e „ (Ci - C )Vu e 



< 



1 



1,1(0) 



\ui. 



lx - C )V(u £7i - U) 



dx 



\ui c 



\-C )VU 



dx 



< CU\\ H -U2 



(an): 



hence, possibly extracting a subsequence, that we do not relabel, we may assume 
that, for some matrix- valued measure rj, 



1 



-Xu e (Ci - <Co)Vu Cn dx ->• dr] 



in the weak* topology of (C°(fi))'. 

We must now show that, for any scalar function $, 



$dr]= / $MVt/d^. 



(3.16) 



In order to do this, it is enough to prove that 



(Co-C^W: Vv i i$dx= / (C - Ci)Vu £ : Vv ij <f> dx + o(|w e 



(3.17) 



because then, by passing to the limit along subsequences of u> Cn in (|3.17|) . we get 
(|3l6]l . 

Let us notice that, since 

( div(c V[/) = div(<C e Vu £ ) in H 
\ (C VC/)t/ = (C £ Vu £ )j/ on dSl, 

for every vector valued test function \I/ we have that 



(3.18) 

For the same reason 
(3.19) 



iVC/ : Wdx 



; Vji £ : V* dx. 



C Vv ij : Wcfe= / C £ V^ J ' : Wdx 
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We can calculate 



f (C - C £ )W : Vvf® dx - /(C - C £ )Vu £ : Vv ij $ dx 
Jn Jn 

: W : V(«* J #) - C e W* : V(U$j) dx 



C W J : V(u £ $) - C £ Vu £ : V(w y $)J 
C W : (vl j ® V$) - C e V< j : (U ® V$)) 



+ y (C W J : (u 6 ® V$) - C £ Vu £ : (v %3 ® V$) ) da: 



By (|3.18p and (|3.19p and recalling that C £ = Co in Q \ uj c , we can write 



(Co - Ci)V£7 : WJ$dx - / (C - Ci)Vw £ : Vv ij ®dx 



= J (C £ Vu £ : V« J $) - C Vu y : V(C/$) ) dx 



C 



C e Vi>^ : V(w £ $) - C VC7 : V(u y $)J dec 
C W : (uf (8 V$) - C e Vu« : (17 ® V$)) dx 
C W J : (u e ® V$) - C £ Vu £ : (w ij ' ® V$)) dx 



£ Vu e : Vi#$ - C W j : W$") dx 



C e Vu^ : Vu £ $ - C VJ7 : Vv lJ $J dx 
C £ Vu £ : (v^ <g) V$) - C W y ' : (U <g) V$)) dx 
C £ V^ J ' : (u £ ® V$) - C VC/ : (i/ J ' ® V$)) dx 
C W : («** ® V$) - C 6 Vv« : (17 ® V$)) dx 
C W J : (w e ® V$) - C £ Vu £ : (w ij ' ® V$)) dx 
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By symmetry of the elasticity tensors, the first two lines of last equality give zero. 
By rearranging the various integral in a suitable way we get, 

f (C - Ci)W : Vvi j $dx - f (C - Ci)Vu £ : S7v ij $dx 
(Vw e - w) : ((v? - v ij ) ® V<f>)dx 
Z c (Vv ij - V^ J ') : ((U - M e ) <E) V$)cfe 



+ / (Ci - C ) VU : ((«« - v lJ ) ® V$)dx 
(Ci - Co) Vv lJ : ((U - u e ) ® V$)dx 
By Lemma 13.11 and by regularity of functions U and in Kq , we get ()3.17[) . 



□ 



3.3. End of the proof of Theorem 12.11 Let Lu Cn as above. By the definition of 
the Neumann matrix it is easy to see that 

(3.20) (u en - U)(y) = f (d - C )Vu e „ : VN(;y)dx. 

(See [4] for details). 

Let Kq C ft the compact set introduced in (|2.8|) . Given y £ 9f2 it is possible to find 

£ C°(fl) a matrix valued function such that ^ y (x) = V x N(x, y) for x £ K . 
Using the previous lemma we get 

{u €n -U)(y) = | We J [ ixo, e „(Ci-C )Vu E „ :*„da: 
Jn l w «J 

= |w e J / MW : %dn + o(\w e \) 
Jn 

= \ui €n \ I MVU :VN(;y)dn + o(\ui e \). 
Jn 

4. Properties of the elastic moment tensor 
In this section we prove few basic properties of the elastic moment tensor M: 

Proposition 4.1. The elastic moment tensorWl has the same symmetry properties 
of the elasticity tensors Co and Ci , that is 

Mtjki = M Hi: j = Mjiki, fi-a.e. 

for any choice of indices i,j, fc, / between 1 and d. 
Moreover, for any symmetric matrix E, 

(4.1) CqC^ x {Cx-Cq)E:E<ME:E<{C x -£, q )E:E \i-a.e.. 

Proof. Firstly, we show that M enjoys the same symmetry as the elastic tensors Ci 
and Cq. To this end, we recall the following equality, which was obtained in the 
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proof of lemma 13.21 

j^-r I (Ci-Co)W: Vv{ j $dx= j^-r I (Ci-C )Vm £ : Vw y $dx + o(l). 

Substituting U and u e for v hk and v^ k respectively, we see that 

-^-r [ (Ci - CojVv" : Vt##da: = ( (C x - C )V^ lfc : V« y *da: + o(l). 
Recalling (|3.14[) and by the symmetry of Co and Ci, we get on one hand that 
■r^r [ (C x - C )Vv hk Vv?*dx = A / y2(Ci-C ) lmpq S hp 6 kq ^^^dx 



Impq 



d(v 13 ) 

1 - C )imfefc V — ®dx 

oxi 



On the other hand 



in 

It follows that 

(4.2) Mijhfe =M hfei j, /i-a.e. 

To obtain the minor symmetry, we observe that 

gf v v ) ) 
2J IXw c „ (Ci — Co); mpg ^* g = |<^eJXw e „ (Ci — Co) m ; pg ^" 9 , 

since Ci — Co is a totally symmetric 4th-order tensor. Then, from p,15|) . it follows 
that 

M i jim=M i j m i, fi-a.e. 
All the minor symmetries now follow from (|4.2p . 

For the proof of (|4.1[) . we follow [6], where the case of the scalar conductivity 
equation was discussed. 

We begin by hxing a constant, symmetric matrix E = [Eij]. We set 



(4.3) V = Y J E l] v L \ v e = J2 E ij 



where v 1 ^ and v l J are given in (|3.12[) and (|3.13l) . and observe that V solves 

Jdiv(C W) = div(Co-B), in ft, 
\ (C VV> = (C E)v, on dfl, 



(4.4) 
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and, consequently, v e solves 



(4.5) 



fdiv(C £ Vv e ) = div(Co^), in O, 
\ (C e Vv e )^ = (CqE)v, on dil. 



We next recall that M is obtained as the weak limit (|3.15[) . Therefore, given any 
function <5 £ C°(f2), we have 



1 



Ja H 



(Ci-Co)Vu^ §{x)dx + 0(1), 
fci 



where e is an element of the sequence {e n }, and square brackets indicate compo- 
nents. We choose $(x) = Eij Eki4>Xni where is a positive smooth function, and 
sum over repeated indices: 



I ME : E<j)dn = V / E i:j Mkjki E k i4>dx 
n ijkl Ju 



Y, E v f(Ci - Co) ■ V««l BjH^da: + o(l) 
Ne .//,, rrr: L J w 



ijkl 



Recall that Vw y = (e* <E) ej ), hence from (|4.3I) . using that E 1 is symmetric, we have 



W = ^2 E v e * ® e j = £ 5 



so that: 



MB : £0 






1 ^ ^Ojfcipg 



(4.6) 



1 — r / (Ci - C ) Vv e : VV<j)dx+ oil) 



(Ci - C ) W : VV<j>dx 



~ [ (Ci - C ) V(« e - V) : W0 rfx + o(l) 



Let us now notice that 



/ (Ci - Co) V(v t - V) : VV<t>dx = - j C £ V(u e - V) : V(w £ - V)cj)dx + R 
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where, by (|4~4|) . (|4~5j) and by Lemma KT\ 

R := [ (Ci - C ) V(w £ - : VV<j>dx + f C e V(u £ - F) : V{v e -V)(f>dx 
1 > t Jn 

(C e - C ) V(w e - V) : W0dz + / C e V(u e - V) : V(w e - 7)0 da; 
i Jn 

C e Vw e : V(w e - V)4>dx - [ C W : V(v e - F)0dx 



C e V?; e : V((u £ - F)</>) da; - / C W : V((u £ - K></>) dx 

Jn 

C e Vw e : ((w e - V) ® V<pdx + / C W : ((v £ - V) (g) V^da; 

= - / C V(w e - V) : ((v e — V)® \7<f)dx — / CiVu e : ((w e — V)® V(f>dx 
Jn Juj e 

(4-7)- o(M) 

By inserting the above relation in (I4.6[) we get 

[ME:E<j)dii = — / (Ci-C ) VV : VV(f>dx 

(4.8) - / C € V{v e -V) : V(v e - V)cj>dx + o{l). 

Jn 

Since C e is strongly convex, and </> > 0, 

/ C e V{v e - F) : V(> e - V)0dx > 



and, hence, 



ME : E(f> d/z < -r^—r f (Ci - C ) W : VV<j)dx + o(l) 



Jn \Ue 
Since the left-hand side does not depend on e we let e — > and get 

(4.9) f ME : E(j)dfi< [ (Ci - C ) W ; VVcj) d/i - / (Ci — Co) E : E<f>dfj.. 
Jn Jn Jn 

Now, by (|4.7I) and by the fact that Ci is strongly convex, we have that 
C e V(w e - V) : V(v e - V)(j)dx 

(C -Cx)V(v e -V) : VV(t)dx + o(\uj e \) 
= I CiV(v e -V) :C^ 1 (C -C 1 )VV<l>dx + o(\u)e\) 
<(^J CxV(v e -V) : V(v e -V)(f>dx\ 

(4.10) (J (C -Ci)W:C7' 1 (Co-C 1 )V^dx^ ' + o(\u £ \) 



1/2 
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from which it follows that 
(4.11) 

V(v e -V) : V(v e -V)c/>dx < [ Cp^Co-C^VV : (C o -Ci)W0 dx+o(\ui e \). 



n 



By inserting (|4.11|) into (|4.8j) we get 

I ME : E(j)dfi > -5-r / (Ci - C ) W : W^tfat 



-r^-r / Ci^Cq-CiJW : (C o -Ci)VV r 0dx + o(l) 

Pel iu E 

= 7^7 / Cr x (Ci - C )W : CoWdz + o(l). 

We | ■/«„. 



By letting e->0we get 
(4.12) 

/ ME:Ecj)dfi> [ CoC^ 1 (Ci-Cq)VV :VV4>d(i= [ CoCr^Ci-Co)^ : E^dfj,. 
Jn Jn Jn 

Notice that (|43jl and (|4.12|) hold for every positive </>, hence (|4.ip follows. □ 



5. The case of thin planar inclusions 

In this section, we specialize to the case of thin inclusions in a planar domain 
SlcK 2 , modeled as an appropriate neighborhood tu e of a given simple curve oo C O, 
that is: 

(5.1) UJ e = {x G il : d(x, (To) < e} ■ 

We impose the following conditions on ctq. We assume that (To is of class C 3 and 
that there exists some K > such that 

d(o- ,an) > r- 1 

(5.2) ||(t || C 3 < X 

R- 1 < lcngth((T ) < K 

Moreover we assume that for every x G (To there exists two discs i?i and Bi of 
radius If -1 , such that 

B 1 n b 2 = Bi n cr = s 2 n <r = {x}. 

The latter assumption guarantees that different parts of cto do not get too close, so 
that ui e does not self- intersect for small e. We refer to (To as the support of uj t . 
Let us fix an orthonormal system (n, r) on a such that n is a unit normal vector 
field to the curve and r is a unit tangent vector field. If a is a closed curve, then 
we take n to point in the outward direction of the domain it encloses. 
We present a different derivation of the small volume asymptotic formula (|2.12[) 
for the displacement at the boundary in this case, which makes more explicit the 
measure and elastic moment tensor M that appear in (|2.12j) . 

The main result of this section is the following theorem, which is a counterpart to 
Theorem [2~T1 
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Theorem 5.1. Let fl C R 2 be a bounded smooth domain and let ao CC be a 
simple curve satisfying JJj|). Let u e and U be the solutions to \2.1$ and [2.7\ l 
respectively. For every x G ao, there exists a fourth order elastic tensor field M(x) 
such that, for y £ <9f2 

(5.3) (u e - U)(y) = 2e M(x)Vu {x) : VN(x, y) da (x) + o(e). 

J (TO 

The term o(e) is bounded by Ce 1+ \\ip\\ H -i/2^gQ^, for some < 8 < 1 and C de- 
pending only on 9, fl, ao, (3o and K . 

5.1. Proof of Theorem 15. li The proof of the theorem closely follows the proof 
of the corresponding result in the isotropic case (see [4]). We only detail those 
steps, where the proof differs from that case. In the following we set u\ = u e ^ 
and u\ = u e | s1 ^ ■ We simply use u e when no confusion can occur. Firstly, we 

write (ii e — U)\ gn in terms of an integral over ui e of the product of Vu* and ViV. 
Secondly, using some regularity estimates for solutions to the elastic system in a 
laminar domain due to Li and Nirenberg |10) . we approximate this integral by 
an integral over a portion of dui t , which we rewrite, in a third step, using the 
transmission conditions and a tensor M satisfying (I5.7[) . The existence of M is 
proved later in Subsection 15.21 Finally, taking limits in the resulting expression as 
e — > and using fine regularity estimates for u e proves the theorem. 
First step. 

We recall (see (j3~20| ) that, for y G dVl 

(5.4) (u e - U)(y) = f (C 1 -C )Vul:VN(-,y)dx. 

Second step. Let f3 be a constant, < j3 < 1, and set 

co' £ — {x + t n(x) : x G <tq, d(x, dao) > e' 3 , t G (— e, e)} . 
Notice that if ao is a closed simple curve, then uj'^ = co e . 

By Theorem 2.1, chapter 2 in combined with Sobolev Embedding Theorem, we 
have that ||VC/||£°o( We ) and || ViV(-, 2/)||l=( W£ ) (for y £ <9f2) are bounded uniformly 
in e. Using this fact together with the energy estimate (|3.3p . one can easily show 
as in [3] that 
(5.5) 

(Ci - Co) Vul : VN{; y)dx= [ (Ci - C ) Vu' : VN{; y) dx + 0(e 1+/3/2 ). 



Let a' e denote the curve 

a' t — {x + fn(i) : x € ao, d(x, dao) > e' 3 } . 

A crucial ingredient, at this point, is a C a regularity estimates for the gradient 
of solutions to laminated systems due to Li and Nirenberg (see [10]). Using this 
estimate and proceeding as in [3], we can approximate the values of Vu £ in uj e ' by 
its values on a' e , so that 



(Ci - C ) Vu l e {x) : VN(x, y) dx 
(5.6) =2e [ (d - Co) V u\ : X7N(-,y) + 0(e 1+a - 1 3 {1+a} ), 
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for /3 < a(l + a) -1 . 
Third step. 

From the results of the section 15721 for every i £ cr|, there exists a fourth-order, 
symmetric tensor M(cc), independent of e, such that 

(5.7) (d(x) - Co(a:)) Vu*(x) = M(a:)Vi£(a:) 

Inserting ()5.7|) into (|5.6|) . we get 

(5.8) 

(C 1 -C )Vu i e :VN(-,y)dx = 2e / MV< : ViV(-, y)da: + 0(e 1 - a -^ 1+a >). 



Fourth Step. 

Now we show that 

(5.9) UV^-WIUoo^CelVII 
for some positive 7. 

Once estimate (|5.9I) is proved, then (|5.8p holds with Vf7 instead of Vuf and (I5~ 
follows immediately by continuity. 

In [4] the proof of estimate (I5.9P strongly relies on the special features of a homo- 
geneous and isotropic tensor Co- In the present case, we use of a Caccioppoli-type 
inequality proved in the Appendix. 

Let 2e < d < d /2 and Cl e d = {x £ Q : d{x, d{Q\u; e ) > d}. 
Since u e — uq is solution to 



div ^CoV(u £ — uq)J =0 in Q \ u) e , 



the regularity assumption on Co implies that u e — uq € iJ ; 2 oc (il\u; e ) ( see [5j Theorem 
2.1 chapter 2]). 

Let k € {1, 2} and let </> £ := dk{u t — uq). The function <f> k solves 

div(CoV^) =F in fi\w e . 

with 

F = -div ((a fe C )V(u £ - uo)) • 
By Caccioppoli inequality (Theorem I A. 1 1 for u = 0) and by p.3[) , we see that 

llVtfll^nj,,) < §H« 2( n* /4 ) + \\F\\l-^ d/d 

= C (^H^Hl^J + IIVK - U o)|li 2( a W) ) 
< CW\%-i/z (d~ 2 + 1) e < C||^||^- 1/2 ^ 2 e. 

and, by (|3.3|) again 

ll^llffi((n« /2 )) < C||V||H-i/,d- 1 Vi 
Theorem IB . 1 1 applied to fc shows that 

||V0 fe || i2+ ,(^) < C (iIFUh-x,^^^) + d^-i||V^*|| i3(n . /a) ) . 

Notice that, by applying Theorem IB. II to « e — u in 17 \ cj £ , 

H-FlU-w^n^) = ||(<9/cC )V(w e - uo)||x,a+n(n« /2 ) 

< C||V(u 6 -uo)||z,2 + r, (n = /2 ) < ||V(« £ - uo)IU=(n« /4 ) 
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and, hence, 

||V/|| i2+ „(^) < c(||VK-n )|U2(^ /4 )+^- 1 ||V/|U 2( ^ /2 )) 

< C (d^^ 2 ) y/l 



On the other hand, applying Theorem IB. II to u t — uq shows that 



By Sobolev Embedding Theorem, it follows that, for k = 1,2, 

(5.10) \\d k (u e - uo)IU-(n 5 ) = \\<t>eh°°(ni) < C (d^- 2 ) ^ 

Now, let y € da' e and let yd denote the closest point to y in the set £l d . From the 
gradient estimates for u e and uq (see [10] and [H Prop. 3. 3]), we have 

(5.11) |V<0/) - V<(y d )| < Cd a , 
which yields, by (|5.10[) . 

|V« - uo)(y)| < C(d a + d-z+dte 1 / 2 ) . 



Choosing d = e 2(2+a 5 we get 

(5.12) |V«-n )(y)|<Ce^ 

where 7 = 2 — 5— r- and hence 

2 ( 2 +"-2+^) 

\\V{ul-uv)\\ L ~ (<) <C^ . 

We conclude exactly as in [4] by noticing that the tensor M is continuous in fl (see 
next section). □ 



Remark 5.2. If the elasticity tensor Co is smoother than C 1,a in ft, by differenti- 
ating again the equation for u e — uo in Sl\w £ and using again Caccioppoli inequality 
we obtain a better exponent"/ in i5.12\) . Moreover, if Co £ C°° this result can be ex- 
tended to any dimension for small neighborhoods of regular hypersurfaces contained 
in SI. 

5.2. Construction of the elastic moment tensor. To establish the asymptotic 
formula (|5.3j) . one seeks to express the term (Co — Ci)Vu* (x) as a linear function of 
Vit'(i), for x <E cr' e (with the notation introduced in the second step of the previous 
proof). This linear dependence defines the elastic moment M in this context, which 
therefore must satisfy: 

(5.13) (C - Ci)W e (x) = M(x)V<(x) 

Let us also recall that n(x) and t(x) denote the normal and tangential directions 
to a' e at point x. 

From the discussion in the previous subsection, it suffices to obtain M along a' e , 
away from the vertices of the curve (more precisely, outside a disk of radius e 13 
centered at the vertices). 

We introduce a coordinate system in ui' e , adapted to the geometry of the problem, 
as follows. Let r be a unit tangent vector field along the curve a' that supports u t , 
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for instance the velocity vector field using arclength, and let n be a vector field 
normal to ctq so that {r, n} forms a frame field on cto. 

We employ again the frame field (r, n) along a' , and extend it to u>' e by setting 
t(x') and n(x') to be constant along the segment {x = x' + hn(x'), < h < e}. By 
construction a global coordinate system in io' e is given by x = x' + hn(x'), where x' 
is a point along the curve o~' and < h < e. We can extend this coordinate system 
smoothly to part of the boundary, that is, for h = e. 

Since the part of the boundary tu e above o~' is given by the graphs of two smooth 
functions above the same curve, by the regularity results recalled in the previous 
section (see [10]), the interior and exterior strain fields are regular and satisfy the 
following transmission conditions pointwise on a' e : 

(5.14a) V<(:e)t(x) = V?4(x)t(x) 

(5.14b) C Vu e e (x)n(x) =CiVu{(i)n(i). 

A direct determination of a 4th-order tensor M that satisfies (|5.13|) directly from 
the above relations turns out to be rather complicated. One is led to invert a linear 
system, the determinant of which is not easily seen to be non zero when Co and Ci 
are anisotropic tensors, unless additional assumptions are made. 
Instead, we follow a construction due to G. Francfort and F. Murat [8]. In their 
work, effective elastic properties of laminate composites are calculated in terms of 
the geometric parameters of the phase layers (see also [TT], P. 168). This construc- 
tion precisely relies on the transmission conditions (|5.14[) . We fix a point x' G o~' , 
and denote the values of the interior and exterior strains respectively by 

e' = V<(i) and e e = V<(x). 
Proposition 5.3. There exists a 4~th order tensor M such that 

(5.15) (Co-Ci)e* = Me e , 

where M depends on Co, Ci, n, and r, but not on e e and e l . 

Proof. First of all let us notice that, by the first transmission condition (|5. 14aft we 
can write 

(5.16) e e = e % + S <£> n + n <g> 5, for some 6 e R 2 

Let q~ x denote the (symmetric) linear operator associated to the quadratic form 
on K 2 : 

{q~\)i = (d(C®n)) : (€®n), (,^t 2 . 
Note that, since Ci is positive definite, q~ x is invertible since, by (|2.5p , 

V(£R 2 , (q-'OC = (Ci(C®n):(C(8n) > A |C| 2 , 
and thus its inverse q is well defined. Conditions (|5.14b[) and (|5.16[) imply that 

(5.17) ((Co -Ci)e e )n) = Ci(e l - e e )n 

(5.18) = (Ci(6®n + n®5))n 
On the other hand, for any £ £ K 2 , we have 

S^^K = 2(Ci(tf®n)) : (£®n) 

= (Ci (S ® n + n ® 5)) : (£ ® n) 
= Ci(e l -e e ) : (£(g)n). 
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This can equivalently be written as 

2(q- 1 6)£= (d(e 4 - e » V£eK 2 

that is, by ([5~T7| 

2(?- 1 <J) = d(e l - e e )n = ((C - C x ) e e )n. 
From the above equation we deduce that 



5=ig((C -C 1 )e e n) 



hence, by (|5.16p 



e l = e e + ~g ((C - Ci)e e n) ® n + n <g> ~g ((C - Ci)e e n) . 

We can now conclude that 

(Co - d) e * = (Co - d)e e + (Co - d) (g ((C - d)e e n) <g> n) , 

where, in the last expression we used symmetry of tensor Co — d- 
Hence, the forth order tensor M defined by 

(5.19) Mh = (C - Cx)h + (C - Cj) (g ((C - C^hn) ® n) , 

satisfies (|5 . 15|) . As can be seen by its expression, M does not depend on e e and 
e 1 but only on the elasticity tensors Co and Ci and on the directions n and t. 
Moreover, this tensor M can be defined for every point x in ujJ and it is continuous 
with respect to x. 

Remark 5.4. Assume that {uj e } is a sequence of thin strip-like inclusions as 
in \5.1)) . Then, the expansions i2.12\) and &5.3\) coincide. In other words, the 
measure M.(x)dpL x that appears in H2.12\) is precisely M.(x)6 ao (x), where M(x) is 
defined in the section \5.2\ 

Indeed, given the form of the sets U) e , it is immediate to see that 
Thus, recalling Lemma \3.2l it suffices to show that 

r^TXaJCl-CoJVUe, M(x)VU(x)5„ , 

in the weak* topology of [C (ft))'. This again is a consequence of the uniform 
regularity estimates on u e . Indeed, let <fi G C (Cl). Using the same notations and 
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the same analysis as in section [5J\ we see that 



-(Ci — Co)Vw e 0dx 



-(Ci - C a )Vu e (j)dx 



■(Ci — Cq)^ u e <pdx 



-> 2 



2^- / (d-CoJVu^do-x + OC^ -^ 1 ^^) 
Vi(as)V«^d<Tx + 0{e int{a ~ l3{1+a) ' l3) ) 
(x)VU(j> d/jt, 



which proves the claim. 



Remark 5.5. In the case of isotropic materials, an explicit formula for M can be 
obtained from (|5.19l) . In fact, in this case Co and Ci have the simple form: 



(5.20) 



-o = A h ® ^2 + 2^4, 
-i = Ai I 2 ® ^2 + 2 mi h 



where Xi, fjb%, i = 0,1, are the Lame parameters, I2, I4 are the identity elements 
on 2 and 4 symmetric tensors respectively. Hence, as linear maps over symmetric 
matrices, they are diagonal in any basis. These tensors are strongly convex if Xi > 
and Xi + fii > 0, i = 0, 1. 

Then, it readily follows that the matrix q is given by (see 8, Equation 4. 8] J: 



(5.21) 



1 



Ai + Mi 



Mi Mi (Ai + 2/ii) 
Combining (|5.20j) with (|5.21|) gives: 

X -X 



n 65 n. 



(g ((C Q - Cx)/in) ® n) = 

2 (/xp - /xi) 
Mi 



— Tr(h) n® n 
Ai + 2mi 

/t \ Ai ~t- Mi //, \ \ 

(/in) (X> 71 (("- n ) • n)n®n 

Ai + 2mi 



Lengthy, but straightforward calculations yield: 

(5.22) M/i = a 2V(/i) 7 2 + &/l + c ((/it) • r) r ® t + cZ ((hn) ■ n) n <g> n, 

w/iere 



a = (A - Ai) 
c= (Mo - Mi) 



A + 2mo 



= (Mo - Mi) — , 
Mi 



Ai + 2mi 

2Ai(mi - Mo) + Mi(Ai - A ) + 2mi(mi - Mo) 



Mi(Ai + 2mi) 



d = 2 (mo - Mi) 



M1A0 — A1M0 



Mi(Ai + 2mi) 

As expected, this formula agrees with that obtained in [4j Theorem 2.1]. 
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Appendix A. A Caccioppoli type inequality 

For the sake of completness, we state and prove a Caccioppoli type inequality for 
solutions of strongly convex systems. 

Theorem A.l. Let B p and B 2p be two concentric balls contained in CI and let u 

be any constant vector. 

Let u <E H (CI) be solution to 

V • (cVu) = f m Cl, 



where C G L°°(d) is a strongly convex tensor and / G H (CI). Then 
(A.l) || Vu\\l HBp) <^\\u- u\\h (B2p) + C 2 ||/||!_ 1(fl2p) 

Remark A. 2. Note that u = (u)b 2p := i Bs i J B2 udx gives the minimum value 
for \\u~u\\ 2 L2{B2p) . 

Proof of Theorem \A.1\ 
Let G C^(B 2p ) such that 

C 

< < 1, = 1 in B p and |V0| < — . 

Let v G H 1 (B 2p ,R 2 ) and observe that 

V(0 2 u) = 9V(9v) + 0V0 ® v. 
By assumption u is a weak solution of 

V ■ (cVuj = f in B 2p 

hence we have 

a Q (u,4>) = -<f,4>> V0eH'(B 2p ) 

where 



o o (u,0) = / CVu : V</>. 
Let ■0 = 9( u — u) and let 4> = 9i\). We have that 
a (^) = / CW>: ViP 



[C(0V(u - u) + (u - u) ® V0) : (0V(u - <S) + (u - u) <g> V0)] 
After some straightforward calculations one gets 

a (ij), V) = a (u, <j))+ C (V0 ® (u - u)) : (V0 ® (w - tt)) =<f,<p> +1 



We have 

|/| < C* 3 ||V0 ® (« - S)|||2 (B2p) < ^||« - u||i» (fla(() 



and by Young inequality 

!</>>!< II/IIh- 1(B2p )I|V0|| L 2 (B2p) < <5||V0|| 2 L 2 (B2p) + c (^)||/|| 2 ff - 1(B2p) 
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So, we get 

a o (^V)-<5||V0||| 2(B2p) < ^\\u-u\\l 2{B2p) +c(5)\\f\\ 2 H - 1{B2p) 

By the strong convexity of C and by Korn inequality applied to the test function 
t/j, we have 

ao(i>,ip) > 7llW|li=( B2p ) 

which gives 

^t\\u ~ u|||a (B2p) + c(S)\\ff H - HB2p) > J\\ V^||| 2(B2p) - 5||V<^||| 2(B2p) 

> j\\^\\h {B2p) -s\\vm\\i HB2p) 

Hence, since 6 < 1 and |V0| < — , 
C 

^■||«-«lli»(B ap )+c(*)||/|||r-i(B ap ) 

> 7l|V^||| 2(B2p) - 28\\Vi>\\ 2 LHB2p) - fWWl'fr,) 

> 7l|V^|| 2 L2(B2p) - 25||V^||| 2(B2p) - -j||u - u||| 2(B2p) 

Choosing 5 small enough, using the fact that 9 < 1, that | V#| < — and the definition 
of "0 we finally get 

l|Vu|| 2 L2(Sp) = ||V^||l a(Bp) < ||V^||i 2(B2p) <^\\U- u\\h( B2p) + C2\\f\\ 2 H - HB2p y 

□ 



Appendix B. A Meyer's type result 

We state a generalization of Meyer's theorem concerning the regularity of solutions 
to systems with bounded coefficients. For r\ > 0, define _ff 1,2+,? (0) by 



H 1,2+T? (0) := |u e L 2+,1 (n), Wu E L 2+ri {n) 
and let H^ 1 ' 2+r '(^l) be its dual. Introduce 

h^(si) ■.= lu e H^+^^.yx cc n|. 

Theorem B.l. There exists r\ > smc/i £/ia£ if u G is solution to 



where C G L°°(f2) is a strongly convex tensor and f <E H 1,2+v (fl) then u £ 
Hj' 2 c +n (Cl) and given B p and Bi p concentric balls contained in Q, 

(B.l) ||Vu|| L2+ „ (Bp) < C(\\f\\ H -i,2 +HB2p) + p^- 1 ||Vu|| L 2 (B2p) ). 

Theorem IB . 1 1 has been proved by Campanato in [5] in the case of strongly elliptic 
systems. From the proof of this result in }5, Chapter II, section 10], it is clear that 
the result can be extended to more general systems provided the Caccioppoli type 
inequality fTheorem lA.il) holds. 
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